Exact Casimir interaction between eccentric cylinders 
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The Casimir force is the ultimate background in ongoing searches of extra-gravitational forces in 
the micrometer range. Eccentric cylinders offer favorable experimental conditions for such measure- 
ments as spurious gravitational and electrostatic effects can be minimized. Here we report on the 
evaluation of the exact Casimir interaction between perfectly conducting eccentric cylinders using 
a mode summation technique, and study different limiting cases of relevance for Casimir force mea- 
surements, with potential implications for the understanding of mechanical properties of nanotubes. 
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As the size of physical systems is scaled down into 
the micrometer and submicrometer scales, macroscopic 
quantum effects become increasingly important. Forces 
attributable to the reshaping of quantum vacuum fluc- 
tuations under changes in geometrical boundary condi- 
tions, predicted almost sixty years ago by Casimir []J, 
have been measured in recent years with increasing ac- 
curacy 0; 0, 0, IE IE ■ The fact that the magnitude and 
sign of the Casimir force depend both on the geometry 
and material structure of the boundaries paves the way 
to several opportunities and challenges for engineering 
mechanical structures above the nanoscale 0- 

The Casimir interaction for perfect metals has been ex- 
actly evaluated only for a limited number of geometries, 
starting from the original parallel plate configuration pj . 
Until recently, for all non-planar geometries the Casimir 
force has been estimated using the so-called proximity- 
force approximation (PEA) |j| , semiclassical 0, and 
optical approximations [12| , and numerical path-integral 
methods [l3| . This has originated a debate on the assess- 
ment of the accuracy of the measurements, a crucial issue 
to establish reliable limits on extra-gravitational forces 
in the micrometer range, fn the last months large vio- 
lations of PFA for corrugated plates have been reported 
|l4j , and the exact Casimir interaction between a sphere 
in front of a plane and a cylinder in front of a plane has 
been computed 0,0,0]. As first discussed in 0|, the 
cylinder-plane configuration is intermediate between the 
plane-plane and sphere-plane geometries, offering easier 
parallelization than the former one, and a larger abso- 
lute signal than the latter one due to its extensivity in 
the length of the cylinder. A related experimental at- 
tempt aiming to measure temperature corrections to the 
Casimir force is under development |19| . 

In this Rapid Communication we present the exact 
evaluation of the Casimir interaction for another geome- 
try of experimental relevance consisting of two perfectly 
conducting eccentric cylinders. Although parallelism is 
as difficult as for the plane-plane geometry, this geome- 
try offers several experimental advantages. First, Gauss 



law dictates that the expected gravitational force is zero 
for any location of the inner cylinder, which allows for a 
null experiment when looking for intrinsically short-range 
extra-gravitational forces. Second, the fact that the con- 
centric configuration is an unstable equilibrium position 
|l8l | opens the possibility of measuring the derivative of 
the force using closed-loop experiments. Finally, residual 
electrostatic charges on the surfaces can be exploited to 
maximize the parallelism between the cylinders looking 
at the minimum value of the resulting Coulomb force. 

Before embarking on the exact calculation of the 
Casimir energy for this geometry, let us recall the result 
of the PFA. This is a simple, though uncontrolled way, of 
treating non-planar configurations, and is valid for sur- 
faces whose separation is much smaller than typical local 
curvatures. For two very long eccentric cylinders of radii 
a < b, length L, and eccentricity e (see Fig. 1), the 
PFA approximation for the non-concentric Casimir en- 
ergy reads £pfa = — 7r 3 fi,cLe 2 /120a 4 (o: — l) 5 , valid when 
a = b/a — > 1 and for small eccentricity S = e/a <C 1. In 
the limit of large eccentricity (e sa b — a) PFA predicts a 
behavior similar to that of a cylinder in front of a plane 
|l8l |. The geometrical dimensionless parameters a and S 
fully characterize the eccentric cylinders configuration. 

In order to go beyond the PFA result, we start by ex- 
pressing the Casimir energy as E = (^/2)X! p ( w p — 
where lo p are the eigenfrcqucncics of the electromag- 
netic field satisfying perfect conductor boundary condi- 
tions on the cylindrical surfaces, and u> p are the corre- 
sponding ones to the reference vacuum (cylinders at in- 
finite separation). In cylindrical coordinates, the eigen- 
modes are h n ^_ = i?n(r, 0) exp[— i(u> n ,k z t — k z z)], where 
uj p = uj n- k z = \Jk"l + A 2 , and R n (A„) arc the eigen- 
functions (eigenvalues) of the 2D Hclmholtz equation. 
Using the argument theorem the sum over eigenmodes 
can be written as an integral over the complex plane, 
with an exponential cutoff for regularization. In order 
to determine the part of the energy that depends on the 
separation between the two cylinders it is convenient to 
subtract the sclf-cncrgies of the two isolated cylinders, 
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FIG. 1: Geometrical configuration studied in this paper. Two 
perfectly conducting eccentric cylinders of radii a < b, length 
L, and eccentricity e interact via the Casimir force. The equi- 
librium position at e = is unstable - any small perturbation 
will make the cylinders snap into each other. 



E 12 (a, b,e) = E — E\{a) - Ei(b). Then the divergencies 
in E are cancelled out by those ones in Ei{a) and Ei(b), 
and the final result for the interaction energy is 



, , , hcL r 

E 12 (a,b,e) = —J o 



dy y log M(iy). 



(1) 



Here M = ( J F/F 00 )/[(F 1 (oo)/F 1 (a)) (F 1 (oo)/F 1 (b))}. 
The function F is analytic and it vanishes at all the eigen- 
values A n (i^oo, at A„), and, similarly, F± vanishes for all 
eigenvalues for the isolated cylinders. The function M is 
the ratio between a function corresponding to the actual 
geometrical configuration and one with the conducting 
cylinders far away from each other. As this last config- 
uration is not univocally defined, we use this freedom to 
choose a particular one that simplifies the calculation. It 
is convenient to subtract a configuration of two cylinders 
with very large and very different radii, while keeping the 
same eccentricity of the original configuration. Eq. Q 
is valid for two perfect conductors of any shape, as long 
as there is translational invariance along the z axis. 

The solution of the Helmholtz equation in the annu- 
lus region between eccentric cylinders has been consid- 
ered in the framework of classical electrodynamics and 
fluid dynamics pol l2lj. The eigenfrequencies for Dirich- 
let boundary conditions (TM modes) and for Neumann 
boundary conditions (TE modes) are given by the zeros 
of the determinants of the non-diagonal matrices 



Q 



TM 

mi i 
TE 
mn 



[J n (Xa)N m (Xb) 
[J' n (Xa)N^(Xb) 



J m (Xb)N„(Xa)} J„_ m (Ae), 
J' m (Xb)N' n (Xa)} J n _ m (Ae), 



where J n and N n are Bessel functions of the first kind. 
The function M can be written as M = M TE M™, 
where M is built with (R being a very large radius) 

F™ = det [Q™(a, b, e)Q™(b, R, 0)] J[ J n (Xa), 

n 

F™(a) = det [Q™( 0) R, 0)] FJ J n (Xa), (2) 
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FIG. 2: Ratio of the exact and PFA Casimir interaction en- 
ergy differences AE — E\i — between eccentric (E12) 
and concentric (EI2) cylinders in the limit of small eccentric- 
ity e <g a. The curve EM denotes the full electromagnetic 
Casimir energy. 



Similar expressions hold for M TE . 

The Casimir energy can be decomposed as a sum of 
TE and TM contributions 
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with M TE <™(^) = det[6 np 
diagonal matrices Aj® and A™ are 



/(TEjTMl 

np J 



(3) 

The non- 
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np 



A 1 ** - 
rip 



I n (/3) \^ K m (a0) 
K n (f3) ^ I m (a(3) 

j;(/3) v K> m (a(3) 



In 



- n (f35)I m „ p (f3S), 
- n (f35)I m „ p (f3S). 



Here /„ and K n are modified Bessel functions of the first 
kind. The determinants are taken with respect to the 
integer indices n,p = —00, ...,00, and the integer in- 
dex m runs from — 00 to 00. Eq.© is the exact for- 
mula for the interaction Casimir energy between eccentric 
cylinders. This formula coincides with the known result 
for the Casimir energy for concentric cylinders (e = 0). 

„ m , in this particular case the matrices 



is partn 

Bill 



As I n - m (0) 

A™ ,TM become diagonal 

We can also obtain the exact Casimir interaction en- 
ergy for the cylinder-plane configuration [T^ | as a limiting 
case of our exact results for eccentric cylinders Eq.®. 
Indeed, the eccentric cylinder configuration tends to the 
cylinder-plane configuration for large values of both the 
eccentricity e and the radius b of the outer cylinder, keep- 
ing the radius a of the inner cylinder and the distance d 
between the cylinders fixed. Using the addition theorem 
and uniform expansions for Bessel functions it can be 
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FIG. 3: Exact interaction energy difference | AE\ as a function 
of a for different values of <5. Energies are measured in units 
of hcL/Aita 2 . These results interpolate between the (a — 1 )~ 5 
behavior for a — ► 1, and the (a 4 log a) -1 behavior for n>l. 



FIG. 4: Exact interaction energy difference AE as a function 
of 5 for different values of a. Energies are measured in units 
of hcL/A-KQ 2 . The maximum at S — shows the instability of 
the equilibrium point. 



proved that, for 

K m (x + h) 



l„Ax 



(x + h) 

r <' m {x + h) 1 
(x + h) 



l (x)/ p _ m (x) m -K n+p (2h). 



Using these equations (with a; = /3e/a and = (3H/a) 
in our exact formula the known result for the Casimir 
energy in the cylinder-plane configuration is obtained [TtI ] 



-iTM,c-p 
np 



K n+p (2f3H/a), 



A™.<="P = -I^l K n+p (2pH/a). 
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We now focus on quasi-concentric cylinders, when the 
eccentricity is small as compared to the radius of the in- 
ner cylinder, i.e. 8 <C 1. The ratio between the radii, 
a = b/a, need not be close to unity, so that PFA is 
in general not valid in this configuration - as only for 
a — ► 1 do we expect to recover PFA. The zeroth order 
e = corresponds to the concentric case, where the ma- 
trices AF? 1 '™ are diagonal. This case was studied in 

by means of exact and semiclassical treatments, and 
it was shown that the PFA concentric energy goes as 

= — haT 3 L/360a 2 (a — l) 3 . Obviously, symmetry ar- 
guments imply that there is no net Casimir force between 
the cylinders in the concentric configuration. To next or- 
der the correction to the energy, AE = E12 — E^, de- 
pends on e 2 , leading to an unstable equilibrium position 
for the concentric geometry. The behavior of the Bessel 
functions for small eccentricity, I m _ n (z) ~ (z) n_m , sug- 
gests that one should use the tridiagonal version of the 
matrices A™'™ considering only elements with p — n 
and p = nil. Expanding the determinants to order 



0(e 2 ), one can write the TM contribution as 



AE 



TM 



hcLe 2 
' Ana 4 

T>„ - 



n— — 00 



dp f3 3 



1-2} 



n+l 



(4) 



where V™ = I n (P)K n (a(3)/K n (/3)I n (ap) is the e = 0, di- 
agonal TM contribution, and T> n and Af n are the 0(e 2 ) di- 
agonal and non-diagonal TM contributions, respectively. 
They read 



2 4K n (P) 
I n (0)I n+ i(0) 



K n _i(a(i) K n+1 (a/3) 
I n -x(a/3) I n+ i(a/3) 

K n {ap) , K n+1 (aP)~ 



/„(«/?) I„+i(a/3) 



A similar expression holds for the TE contribution, with 
the Bessel functions replaced by their derivatives. Eq.QJ 
and the corresponding TE one are exact expressions for 
the Casimir energy difference between eccentric and con- 
centric cases in the limit e <C a. 

= AE™/2 = 



The PFA limit AE™ A = 
-7r 3 hcLe 2 /240a 4 (a - 1 ^ 5 



A pTM 
PFA 



1)° can be obtained from Eq.QJ 
considering a — ► 1. In this limit the leading contribu- 
tion arises from large values of the summation index n, 
for which the use of asymptotic uniform expansions of 
the Bessel functions is in order. Fig. 2 depicts the ra- 
tio of the exact Casimir energy AE and the PFA limit 
for the almost concentric cylinders configuration. As ev- 
ident from the figure, PFA agrees with the exact result 
at a few percent level only for a very close to unity, and 
then it noticeably departs from the PFA prediction. In 
the opposite limit, a — > 00, the Dirichlet contribution 
is much larger than the Neumann one, and the integral 
and sum in Eq.Q) is dominated by the n — term. The 
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asymptotic result for the energy is 



hcLe 



87ra 4 a 4 log a 
3.3348frcLe 2 
87ra 4 a 4 log a 



d(3/3 3 



K (J3) , K^p) 



HP) h(P) 



(5) 



This equation is valid when log a 3> 1. From Eq.(J5J we 
see that the force between cylinders in the limit a, e <C b 
is proportional to Le/b 4 ln(6/a). The weak logarithmic 
dependence on a is characteristic of the cylindrical ge- 
ometry. A similar weak decay at large distances has been 
found for the cylinder-plane geometry |T^ . 

Next we consider arbitrary values of the eccentricity. 
In this case we need to perform a numerical evaluation 
of the determinants in Eq.JSJ). We find that as a ap- 
proaches smaller values, larger matrices arc needed for 
ensuring convergence. Moreover, for increasing values 
of 6 it is necessary to include more terms in the series 
defining the coefficients ^4™'™. In Fig. 3 we plot the 
interaction energy difference \AE\ as a function of a for 
different values of S. These numerical results interpolate 
between the PFA and the asymptotic behavior for large 
a, beyond the quasi-concentric limit. Fig. 4 shows the 
complementary information, with the Casimir energy as 
a function of S for various values of a, showing explicitly 
the instability of the concentric equilibrium position. 

The exact evaluation of the Casimir force between par- 
allel eccentric cylinders obtained here lends itself to a va- 
riety of applications for experiments implementing this 
geometry in the micrometer and the nanometer scales. 
Firstly, the accurate knowledge of the Casimir force in 
this configuration allows to look for extra-gravitational 
forces while the usual gravitational force, apart from bor- 
der effects, is cancelled. The cancellation of the Newto- 
nian gravitational force is also common to the parallel 
plate geometry if a dynamical measurement technique is 



adopted for the latter configuration, however in addition 
the concentric cylinder case has better shielding from the 
electrostatic force due to spurious charges, allowing for 
a null experiment with respect to both Newtonian and 
Coulombian background forces even with a static mea- 
surement technique. This suggests a micrometer version 
of experiments performed using the concentric cylinders 
configuration and a torsional balance to test the inverse- 
square gravitational law in the cm range p3[ . In particu- 
lar, in the case of a repulsive Yukawian force, one can en- 
visage a situation where the unstable equilibrium due to 
the Casimir force is balanced or overcome by the former, 
and even the qualitative observation of mechanical stabil- 
ity in the concentric configuration will establish the ex- 
istence of a new force. Secondly, the measurement of the 
Casimir force in various configurations involving cylin- 
ders (with the first example provided in £| for crossed 
cylinders) is interesting in itself as one can make reliable 
tests of the PFA approximation versus both the exact so- 
lution and the actual experimental outcome. We expect 
that the case of two parallel cylinders with distance larger 
than the sum of their radii will experience stronger devi- 
ations between PFA and the exact solution. This case 
could be analyzed using the approach presented here. 
Thirdly, our results could lead to a quantitative expla- 
nation, once finite conductivity and temperature effects 
will be taken into account, for the observed easiness to 
bend laterally multiwall nanotubes. The Casimir force 
and its non-retarded counterpart |24| provide a natural 
mechanism to allow for instability from the perfectly con- 
centric situation of nanotubes, and could also play a role 
in the observed decrease of the effective bending modulus 
of nanotubes with large (tens to hundreds nm) radii |25j , 
and their fragmentation [2(j. 
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